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The timing of crucial events in plant life cycles is shifting in response to climate change. We use phenology 
records from PlantWatch Canada 'Citizen Science' networks to study recent rapid shifts of flowering 
phenology and its relationship with climate. The average first flower bloom day of 19 Canadian plant species 
has advanced by about 9 days during 2001-2012. 73% of the rapid and unprecedented first bloom day advances 
are explained by changes in mean annual national temperature, allowing the reconstruction of historic flower 
phenology records starting from 1948. The overall trends show that plant flowering in Canada is advancing by 
about 9 days per °C. This analysis reveals the strongest biological signal yet of climate warming in Canada. This 
finding has broad implications for niche differentiation among coexisting species, competitive interactions 
between species, and the asynchrony between plants and the organisms they interact with. 

Plant phenological events such as bud burst, flowering, and senescence have received increased attention in 
the light of climate warming 1-5 . Even though changes in various plant phenological events have substantial 
consequences for ecosystem functioning, reproductive events such as flowering time largely determine 
population and community dynamics in future generations, affecting evolutionary processes. Flowering time 
also impacts animals for which pollen, nectar, and seeds are important resources, while earlier flowering also 
implies earlier activity in other processes and phenological events (leaf expansion, root growth, and nutrient 
uptake) that are indicators of environmental cues such as temperature, moisture, and photoperiod. The timing of 
spring flowering and leafing as well as terrestrial primary production of perennial plants is largely controlled by 
temperature across Canada 6 " 9 . 'Citizen Science' initiatives have recently been successful at publicizing phenology 
and encouraging public participation in recording. Examples of these are PlantWatch Canada networks 
(www.naturewatch.ca). Among several plant phenological events, first bloom defined as a plant stage at which 
the first flower buds have opened in an observed tree or shrub or in a patch of smaller plants, constitutes the most 
straightforward recording for scientific use. We have analyzed the first bloom day for a set of 19 plant species (see 
Supplementary Table SI online for list of plant species and form) recorded by several observers across Canada 
between years 2001 and 2012 (see Supplementary Table S2 online for list of the measurement locations). We use 
first flower bloom day records from PlantWatch Canada 10 'Citizen Science' networks to study the flowering 
phenology trend and its relationship with climate. 

In the present study, we have selected only the first bloom day records observed by at least five observers, at a 
minimum of five different locations in order to remove observer bias. The interannual variability expressed as an 
anomaly, and trend expressed as a regression slope are calculated for each site and species independently. The 
statistics of interannual variability and trend at species and national levels were calculated from the site-year 
measurement weighted averages of site level anomalies and slopes. Calculating interannual variability and trend 
at each site reduces possible biases due to errors in species naming, geographic position tagging, inter- observer 
bias, and the influences of variability in site productivity and latitudinal gradients on flowering phenology. Our 
approach, i.e., computing site level anomaly and trend statistics rather than directly averaging the first bloom day 
records across Canada, is a novel method for partially removing possible biases when the quality of data is 
unknown in 'Citizen Science' records. 

Results 

Our results suggest that on sample size (site-year measurements) weighted averages, flowering days are advanced 
by 0.732 days per year during 2001-2012 (Figure 1). 18 out of 19 species showed earlier first bloom days 
(Figure 1). 7 (37%) Canadian plant species show significantly earlier first bloom day anomalies (p < 0.05, 
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Figure 1 | The measurement weighted average first bloom days per year for 19 plant species across Canada. The first bloom days per year are calculated 
from linear regression slope (p) in the form of Y = a + |3X. 7 is first bloom day of a given year, a is regression intercept, P is regression slope, and X is a 
calendar year. The numbers (No.) correspond to the number of site-year measurements between years 2001 and 2012. The weighted average of Pearson 
correlation coefficient by the number of measurements per site is also given. The measurements from the Canadian 'Citizen Science' PlantWatch data for 
the given species observed by at least five observers, at a minimum of five different locations, with at least five year measurements between the years 2001 
and 2012, totaling 1784 unique site-year measurements were included in the calculation. The weighted overall average slope by the number of 
measurements is -0.732 corresponding to average advance of flowering days per year. The flowering days advanced over Canada during the years 2001- 
2012 by -8.78 (weighted average by number of measurements). Error bars represent half standard error of the regression slopes. *** represents species 
with significantly earlier first bloom day anomalies (p < 0.05, Student's f-test). 



Student's f-test) (Figure 1). Few species showed larger advances 
including Larch (Larix laricina), Coltsfoot (Tussilago farfara), and 
Mayflower (Epigaea repens). The two arctic and alpine plants, 
Mountain Avens (Dry as integrifolia) and Purple Saxifrage 
(Saxifraga oppositifolia) showed a delay or little change in first bloom 
days. Generally speaking, our analysis indicates that first bloom days 
of temperate and boreal flowering plants are advancing whereas 
arctic and alpine flowering plants show a delay or little change. 

During the sampling period excluding year 2012, Canada's annual 
temperature increased by 0.6°C from 1961-1990 average normal 
(www.ec.gc.ca). The entire sampling period had warmer annual tem- 
peratures than the normal years. Our results confirm that climate 
warming advances the first bloom days of early flowering species 
(Figures 1 and 2). The advances of first bloom days are significantly 
explained by changes in mean annual temperatures (£ 2 = 0.73, p = 
0.0009). Figure 2 shows the statistically significant (R 2 = 0.34, p = 
0.04) Canada-wide first bloom day anomalies during 2001-2012. 
The combined winter (December, January, and February) and spring 
(March, April, and May) temperatures explain 71% of the changes of 
the first bloom days (Figure 2). When the winter and spring tem- 
peratures are used as independent regressors of estimate, they 
explain 75% changes in the first bloom days (Figure 3). Statistical 
analysis between the first bloom day and mean annual national tem- 
perature anomalies indicate that during 2001-2011, plant flowering 
in Canada advanced by about 9 days per °C. In an analysis from 
western Canada spanning 1936-1996 9 , flowering events of Aspen 
(Populus tremuloides) advanced by 8 days in 60 years, a slower rate 
than we observed. 

Discussion 

It is difficult to detect systematic or biological patterns such as phylo- 
geny, life-form, geographical distribution, or pollination mechanism 
in most of our flowering phenological datasets because of the limited 



sampling period which precludes robust analyses. However, it is 
evident that Canadian plants experienced rapid and unprecedented 
advances in flowering time during the past 12 years, largely following 
changes in temperature (Figure 2). Our analysis shows that a power- 
ful climate- warming signal observed throughout the 1990s 2 ' 4,11 " 16 was 
also evident in the 2000s (Figures 2 and 3), the warmest decades on 
record. The long time trend of estimated first bloom days from 
winter and spring national temperatures for years 1948-2012 shows 
a smaller regression slope (-0.087) compared to the observed trend of 
2001-2012 first bloom days which has a larger regression slope (- 
0.669) (Figure 3). This confirms that climate warming is not only 
advancing plant phenological events; it is also advancing the events 
progressively (Figure 3). 

The recent rapid and unprecedented advance of flowering time 
has also been reported by other studies 11 although earlier flowering 
phenology studies show a slower advance of flowering days 31217 . 
Early flowers are more likely to be killed by frosts 18 resulting in a loss 
of the year's seed production in turn making them less able to shift 
their range north quickly enough to keep pace with changing vegeta- 
tion/climate zones. This is important for niche differentiation among 
coexisting species 19 and therefore will alter competitive interactions 
between species disrupting ecosystem structure. Earlier flowering 
trends also confirm how climate change might upset the synchrony 
between plants and the organisms they interact with 20 . 

Phenological networks are experiencing a resurgence of interest in 
Canada (www.naturewatch.ca) and around the world, which will 
greatly enhance our capacity to track the biotic responses to changing 
climate. The participation of 'Citizen Science' for large-scale obser- 
vational studies of plant responses to climate change will enhance our 
understanding of biosphere-atmosphere interactions. There is a 
strong reliance on observational data for studying plant responses 
to changing climate because controlled studies 21 reveal warming 
experiments underpredict plant phenological responses to climate 
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Figure 2 | The relationship between the anomalies of first bloom days and mean national temperatures during 2001-2012 over Canada. Temperature 
anomalies are from 1961-1990 normal. Anomaly of first bloom day is calculated for each site separately for a given species and the arithmetic mean of all 
sites and species are taken as a national average per year. The anomalies of first bloom days are statistically significant against calendar years (R 2 = 0.34, 
p = 0.04). a, The relationship between the anomalies of first bloom days and mean winter national temperatures (R 2 = 0.35, p = 0.04). b, The 
relationship between the anomalies of first bloom days and mean spring national temperatures (R 2 = 0.69, p = 0.0008). c, The relationship between the 
anomalies of first bloom days and mean winter and spring national temperatures (R 2 = 0.71, p = 0.0006). d, The relationship between the anomalies of 
first bloom days and mean annual national temperatures (R 2 = 0.73, p = 0.0009). Excluding the year 2012 where the temperature data is not available, the 
overall trend of first bloom day and annual national temperature anomalies shows plant flowering in Canada advanced by about 9 days per °C. 



change. Our analysis together with 21 affirms responses to climate 
change that are predicted using ecosystem models, or experimental 
data, should be recalibrated using observational data. 

Scientific publication is the most validating dissemination of find- 
ings from data collected by 'Citizen Science' participants. Although 
researchers have only recently begun to explore the integration of 
outcomes of 'Citizen Science' data with changing climate, we recom- 
mend this to be at the forefront of future 'Citizen Science' initiatives. 
Evaluating and disseminating the impacts of 'Citizen Science' pro- 
jects on learning how our ecosystems function may ultimately 
increase the chances of project success while also contributing to 
socio-ecological policies of high priority to society. 

Methods 

Plant Watch Canada is a part of national Nature Watch series of volunteer monitoring 
program. These programs are designed to help scientists collect information about 
the environment on spatial scales beyond the reach of most research budgets. First 
bloom day records were collected from 11 provinces and territories in Canada (see 



Supplementary Figure SI online for distribution of first bloom day records across 
Canada). Records are missing for British Columbia and the Northwestern territories. 
Observers collect first bloom days of common species (Figure 4), georeference their 
data with exact geographic coordinates using web-based maps or a global positioning 
system (GPS), and report their data via email, paper-based report or on-line tools 
(www.naturewatch.ca). Among 13,832 spring plant phenology records across 
Canada, we have selected only first bloom days observed by at least five observers, at a 
minimum of five different locations, with at least five year measurements between the 
years 2001 and 2012, totaling 1,784 unique site-year measurements in order to avoid 
observer bias. First bloom day records from cities, large towns, sloping terrain, and 
shaded sites were excluded from the analysis in order to remove the effects of urban 
heat island, topography, and shade on flowering phenology. The first bloom days per 
year are calculated from linear regression slope ((3) in the form of Y = a + f3X (Y is 
first bloom day of a given year, a is regression intercept, P is regression slope, and X is 
a calendar year) for each species and site separately. The flowering days over Canada 
during the years 2001-2012 are then reported as a weighted average by number of 
site-year measurements. Whereas anomaly of first bloom day is calculated for each 
site separately for a given species and the arithmetic mean of all sites and species were 
taken as a national average per year. 

The annual and seasonal time series of national temperatures from 1948 to 2012 for 
all of Canada were obtained from 330 weather stations synthesised by Environment 
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Figure 3 | Estimated 1948-2012 and observed 2001-2012 first bloom day anomalies. Estimated first bloom days are based on the observed relationship 
between the anomaly of first bloom day (FBD) for 2001-2012 normal, and winter ( T w ) and spring ( T s ) mean national temperature anomalies for 1961- 
1990 normal (FBD = a + piT w + p 2 T s , where a, Pi, and p 2 are regression coefficients). The performance of the estimation is: R 2 = 0.75, p = 0.002, 
and root mean square error of 2.3 days. (Inset) the anomaly of the combined winter and spring temperature from 1961-1990 normal. The fitted 
parameters are: FBD = 3.375792908 - 0.79075471 7 T w - 1.770635276 T s . 




Figure 4 | Major Canadian flowering plant species included in the first bloom phenology study. From top-left to bottom-right: Red Maple, Saskatoon, 
Prairie Crocus, Bearberry, Blue-bead Lily, Bunchberry, Mountain Avens, Mayflower, Wild Strawberry, Larch, Aspen, Rhodora, Labrador tea, Purple 
Saxifrage, Lilac, Dandelion, Starflower, Trillium, and Coltsfoot. Courtesy of Plant Watch Canada. 
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Canada 22 . Since weather stations are not evenly distributed across the country, the 
temperatures are first interpolated using Gandin's Optimal Interpolation method 
based on grids covering the entire country 23 . The seasons are defined as follows: 
winter (December of the previous year, January, February), spring (March, April, 
May), summer (June, July, August), and fall (September, October, November). Both 
the seasonal and annual national temperatures were used to study the relationship 
between climate warming and first bloom days. Finally, the winter and spring mean 
national temperatures were used as independent explanatory variables in linear 
regression to reconstruct first bloom days of the past. 
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